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A systematic look at the Very High and Low/Hard state of GX 339-4: 
Constraining the black hole spin with a new reflection model 
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ABSTRACT 

We present a systematic study of GX 339-4 in both its very high and low hard states from 
simultaneous observations made with XMM-Newton and RXTE in 2002 and 2004. The X-ray 
spectra of both these extreme states exhibit strong reflection signatures, with a broad, skewed 
Fe-Ka line clearly visible above the continuum. Using a newly developed, self-consistent 
reflection model which implicitly includes the blackbody radiation of the disc as well as the 
effect of Comptonisation, blurred with a relativistic line function, we were able to infer the 
spin parameter of GX 339-4 to be 0.935 ± 0.01 (statistical) ±0.01 (systematic) at 90 per 
cent confidence. We find that both states are consistent with an ionised thin accretion disc 
extending to the innermost stable circular orbit around the rapidly spinning black hole. 
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1 INTRODUCTION 

X-ray spectra of Galactic black hole candidates (GBHCs) are an 
important tool in the studies of the inner regions of accretion flow 
around black holes (BHs), providing information on both the ge- 
ometry of the accretion disc and on intrinsic physical parameters 
such as BH mass and spin. 

The spectrum can be explained by the combination of a quasi- 
thermal blackbody component caused by radiatively efficient ac- 
cretion through a disc (Shimura & Takahara 1995; Merloni, Fabian 
& Ross 2000), a power-law component due to inverse Compton 
scattering of the soft thermal disc photons in a cloud of hot elec- 
trons or "corona" (Zdziarski & Gierlinski 2004), and a reflection 
component (Ross & Fabian 1993). The latter arises as hard emis- 
sion from the corona irradiates the cooler disc below and results in 
"reflection signatures" consisting of fluorescent and recombination 
emission lines as well as absorption features. The most prominent 
of these "signatures" is the broad, skewed Fe-Ka line observed in 
a number of GBHCs and active galactic nuclei (AGNs, see Miller 
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^ The term quasi-thermal blackbody applies here to the local emission at 
the surface of the accretion disc. Due to opacity effects it is broader than a 
true blackbody spectrum. 



2007 for a recent review) indicative of reflection from the innermost 
regions of an accretion disc. 

In the inner regions of an accretion disc the iron Ka line shape 
is distorted by various relativistic effects such as gravitational red- 
shift, light-bending, frame-dragging and Doppler shifts, with the 
effects becoming more prominent the closer the line is emitted to 
the event horizon (Fabian et al. 1989, 2000; Laor 1991 ). In the case 
of an accretion disc around a non-spinning Schwarzchild BH, sta- 
ble circular orbits can only extend down to the radius of marginal 
stability, rme— 6.0rgWhere rg= GM/(? . This radius depends on 
the spin parameter a/M, and decreases to ~ l-24rg for a maximally 
rotating (a/M fa 0.998) Kerr BH (Thorne 1974). By making the 
standard assumption that the emission region extends down to the 
radius of marginal stability (i.e. Tm = nns) one can obtain an es- 
timate on the dimensionless spin parameter (Bardeen et al. 1972; 
Reynolds & Fabian 2007). 

GX 339-4 is a dynamically constrained [Mbh > 6.0 Mq; 
Hynes et al. 2003; Munoz-Darias et al. 2008) recurrent black hole 
binary (BHB). Its distance has been estimated at 8 kpc (Zdziarski 
et al. 2004). Observations have been made on multiple occasions 
in various spectral states from the "low/hard" to the "very high 
state" (VHS; for a recent review on the different spectral states 
see e.g. McClintock & Remillard 2006). In both the VHS and the 
low/hard state (LHS) GX 339-4 shows a power-law spectra, with 
photon index F ~ 2.5-2.7 and 1.4-1.5 respectively (Miller et al. 
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2004, 2006, hereafter Ml and M2; for a recent analysis of Suzaku 
observation in the "intermediate" state see Miller et al. 2008), as 
well as the presence of a quasi-thermal disc component, usually de- 
scribed by a multicolour disc blackbody model (MCD; Mitsuda et 
al. 1984). In both cases, the fluorescent Fe Ka features have been 
modelled by the addition of a LAOR relativistic line (Laor 1991) 
plus an ionised disc reflection component (PEXRIV, Magdziarz 
& Zdziarski 1995). In this manner. Miller et al. (2004, 2006) mea- 
sured rin~ 2.0-3.0rg and nn~ 3.0-5.0rg, for the VHS and LHS 
respectively. 

It has long been known that the reflection in a BHB system 
cannot be mimicked simply by adding a blackbody component 
to the reflection spectrum from an otherwise cool disc. Compton 
broadening of the iron Ka line is of greater importance in warm 
accretion discs and should thus modify the spectral behaviour of 
BHB compared to that of AGNs. In this paper we undertake a sys- 
tematic reanalysis of the VHS and LHS of GX 339-4, as reported 
by Miller et al. (2004, 2006). We employ the self-consistent reflec- 
tion model developed by Ross & Fabian (2007), where blackbody 
radiation entering the accretion disc surface layers from below, as 
well as the effect of Comptonisation, is implicitly included in the 
model. 

Our method of measuring the spin of stellar mass BH is com- 
plementary to that of McClintock, Narayan & Shafee (2007). They 
use the soft high state when any power-law emission is minimal 
and fit the quasi-blackbody continuum spectrum. Their method, in 
contrast to ours, requires accurate measurements of the mass and 
distance of the black hole. In the following section, we detail our 
analysis procedure and results. 



2 OBSERVATION AND DATA REDUCTION 

GX 339-4 was observed in its VHS by XMM-Newton for 75.6 ks, 
starting on 2002 September 29 09:06:42 UT (revolution 514) and 
simultaneously by RXTE for 9.6ks starting at 09:12:11:28 UT 
(Ml). LHS observations were made by XMM-Newton during revo- 
lutions 782 and 783, for a total exposure of 280 ks starting on 2004 
March 16 16:23:41 TT and RXTE on 2004 March 17 at 12:03:12 
TT, observation 90118-01-06-00 (M2). For the 2002 observation 
the EPIC-pn camera (Struder et al. 2001) was operated in "burst" 
mode with a "thin" optical blocking filter. For the low/hard obser- 
vation the EPIC-MOSl and EPIC-MOS2 cameras (Turner et al. 
2001) were operated in the standard "full-frame" mode with the 
"medium" EPIC optical blocking filter in place. Starting with the 
unscreened level 1 data files for all the aforementioned observa- 
tions, we followed the reduction procedures mentioned in Ml and 
M2. 

In essence, for the XMM-Newton observation of the VHS the 
Observational Data Files (ODFs) were processed using the latest 
XMM-Newton Science Analysis System v 7.1.0 (SAS), with events 
being extracted in a stripe in RAWX (31.5-40.5) vs RAWY (2.5- 
178.5) space. Bad pixels and events too close to chip edges were ig- 
nored by requiring "FLAG = 0" and "PATTERN ^ 4". The energy 
channels were grouped by a factor of five to create a spectrum. The 
standard canned burst mode response file epnjbu23^dY9jhin.rsp 
was used to fit the spectrum. The total good exposure time selected 
was 2.2 ks. Due to the high source flux, background spectra were 
not extracted. RXTE data for the VHS were reduced using the RXTE 
tools provided in the HEASOFT v 5.2 software package. Standard 
time filtering, including the South Atlantic Anomaly and elevation 
angle restriction ( lOdegs. from the earth's limb) resulted in 
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Figure 1. EPIC-MOSl, revolution 782 spectra for various source extrac- 
tion regions and event selection criterion fitted with a simple powerlaw and 
MCD component(see text). Top: From top to bottom: Source extraction re- 
gion in annulus with inner radius 18" single-quadruple pixel, 18" single 
pixel, 30" single pixel, and 50" single-quadruple pixel. Bottom: Data/model 
ratio for the region around the iron line profile with a powerlaw. The energy 
range 4.0-7.0 kcV was ignored during the fit. It is clear that pile-up is only 
significantly present in the annular with inner radius 18" and both single- 
quadruple pixel (top black). 



a net Proportional Counter Array (PCA) and High-Energy X-Ray 
Timing Experiment (HEXTE) exposures of 9.3 and 3.3 ks, respec- 
tively. To account for residual uncertainties in the calibration of 
PCU-2, we added 0.75 per cent systematic error to all its energy 
channels. The response matrix was made by the task "pcarsp". The 
HEXTE source and background spectra were made using the stan- 
dard recipes. Standard canned response were used for spectral fit- 
ting. 

The XMM-Newton data from the 2004 LHS was reduced us- 
ing 5A5v 7. 1.0. As opposed to the EPIC-pn data, the EPIC-MOS 
cameras in the "full frame" mode are more susceptible to photon 
and pattern pile-up. Pile-up occurs when several photons hit two 
neighbouring (pattern pile-up) or the same (photon pile-up) pixel 
in the CCD before the end of a read-out cycle. If this happen the 
events are counted as one single event having an energy equal to the 
sum of all their energies, thus hardening the spectra. In order to in- 
vestigate the effect of pile-up suffered by the EPIC-MOS cameras, 
we used the SAS tool xmmselect to obtain spectra from annular 
regions of inner radius 18", 25", 30" and 50" and outer radius of 
120" centered on the source. The events were filtered by requiring 
"FLAG=0" and "PATTERN ^ 12" (single-quadruple pixel events) 
as well as "PATTERN==0" (only single-pixel events). Background 
spectra were extracted from a 60" circle near the corner of the cen- 
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tral chip of each MOS camera for both "PATTERN ^ 12" and "PAT- 
TERN==0". Response files for each spectrum were created using 
the tools nnfgen and arfgen. The FTOOL grppha was used to re- 
quire at least 20 COUNTS BIN^^. The spectra of four different ex- 
traction region and event criterion for the EPIC-MOS 1 (revolution 
782) observation are shown in Figure 1, fitted with a simple power- 
law and MCD component modified by absorption in the interstellar 
medium (PHABS model in XSPEC). The various parameter were 
tied between the spectra and a normalisation constant was allowed 
to float between them. It is clear from Fig. 1 that pile-up only sig- 
nificantly affects the spectrum created with the source extraction re- 
gion with inner radius 18" and "PATTERN^ 12". All other spectra 
are consistent with the most conservative extraction region (inner 
radius 50", single-pixel events) at energies between 0.7-10.0 keV. 
The overall shape of the spectrum is, however, not significantly af- 
fected by pile-up in all extraction regions and patterns, as can be 
seen in the lower panel of Fig. 1 . 

In order to maximise signal-to-noise and make use of the best 
calibrated response matrix for the LHS, we use the spectra ex- 
tracted from the annulus with inner radius of 18" and single-pixel 
events throughout the analysis detailed in this work. A net expo- 
sure time of 56 and 59 ks was obtained in revolution 782 for the 
£P/C-MOS 1 and 2 camera respectively and 58 ks for each camera 
in revolution 783. For the RXTE data set, the reduction procedure 
involved the use of the tools provided in the RXTE HEASOFT v 6.0 
software package. We used the "Standard 2 mode" data from PCU- 

2 only. The event files and spectra were screened and the back- 
ground and response files created. Systematic errors of 0.75 per 
cent were added to all PCU-2 energy channels. The HEXTE-A 
cluster was operated in the "standard archive mode". Background- 
subtracted spectrum and associated instrument response files were 
created using standard procedures. The RXTE observations resulted 
in net PCA and HEXTE exposures of 2.2 and 0.8 ks, respectively. 

We restrict our spectral analyses of the XMM-Newton EPIC- 
pn data to the 0.7-9.0 keV band. For a preliminary constraint on 
the blackbody temperature of the LHS we use XMM-Newton EPIC- 
MOS data in the range 0.5-2.0 keV as we expect the temperature 
to be low, however for the rest of the analyses XMM-Newton EPIC- 
MOS is used in the range 0.7-9.0 keV, similarly to the VHS, un- 
less stated otherwise. The PCU-2 spectrum is restricted to the 2.8- 
25.0 keV band with an edge at 4.78 keV (r — 0.1) to account for 
the strong Xe L edge. HEXTE spectrum is analysed between 20.0- 
100.0 keV. A Gaussian line at 2.31 keV is introduced when fitting 
the EPIC-pn spectrum due to the presence of a feature at this en- 
ergy that resembles an emission line. This feature is likely to be 
caused by Au M-shell edges and Si features in the detectors (Ml). 
When fitting the RXTE spectra as well as the four spectra from 
the XMM-Newton low/hard observation, (MOSl and 2 for revolu- 
tion 782 and 783), a joint fit is achieved by allowing a normalisa- 
tion constant to float between the various spectra. All parameters 
in fits involving different instruments were tied. XSPEC v 11.3.2 
(Amaud 1996) was used to analyse all spectra. The quoted errors 
on the derived model parameters correspond to a 90 per cent confi- 
dence level for one parameter of interest (Ax^ = 2.71 criterion), 
unless otherwise stated. 

3 ANALYSIS AND RESULTS 

3.1 Fits to RXTE Data: 2.8-100.0 keV Continuum 

We first analyse the RXTE PCU-2 and HEXTE spectra in order to 
constrain the power-law continuum of the two states. By consid- 



ering the simplest power-law plus MCD model, modified by ab- 
sorption in the interstellar medium (PHABS model in XSPEC) 
with an equivalent neutral hydrogen column density fixed at jVh = 
5.3 X 10^^ cm"^ (Kong et al. 2000), resulted in a poor fit for 
both states, with x^/f = 279.5/74 and 275.0/76 for the VHS 
and LHS respectively. In the case of the LHS the addition of a 
MCD component did not affect the fit. Significant residual fea- 
tures are present in the region around the Fe Ka fluorescence 
line. In order to phenomenologically model a disc reflection line 
we initially added a Gaussian emission line and smeared edge 
components (SMEDGE in XSPEC) to the model. This signifi- 
cantly improved the fit with x^/f = 68.7/68 and 88.9/70 for 
the VHS and LHS respectively. The parameters measured in the 
VHS for this model are Tpl = 2.56toil, -Rpl = 2.7lo 5, kT 
= 0.87j:»-°^ keV, i?MCD = 2200^^°°°,' E^^^,, = e.Ofi ikeV, 
FWHM= 2.6inkeV, 7?gauss = Omtofi, EW = 220 ± 

100 eV, i?smedge = S.Sl'j'^y keV, Tsniodge = 0.21q;2' M^mcdgc = 

2.0l'^ Q keV, (where R is the normalisation for each function). The 
equivalent parameters for the LHS are FpL ~ 1 .48 ±0.01, Rpl = 
0.2041° °°^, Sgauss = 5.9l J-^ keV, FWHM = 0.9t J'^ keV, 

i?gausB = 1.7tli X 10-^ EW = 110 ± 40 eV, i^smedge = 

7.1+'''^ keV, wdgo = 0.3tE!;?, M^medge = 2.3 ± 2.0 keV. The 
values obtained for both states are in agreement with those in Ml 
and M2. 

3.2 Fits to XMM-Newton EPIC-pn and MOS data 

3.2.1 Verifying the presence of a qiuisi-blackbody component 

Whilst fitting XMM-Newton data for both VHS and LHS, the 
power-law index was constrained to lie within AF ^ 0.1 from 
the values obtained in the RXTE fits. We began by considering 
a simple power-law continuum plus blackbody component in the 
form of a MCD. The hydrogen column density was fixed at A'h = 

5.3 X 10^^ cm~^ for both states. A fit to the VHS EPIC-pn data 
in the range 0.7-3.0 keV {x^/u = 1448.9/462) shows the pres- 
ence of a blackbody component with temperature of ~ 0.76 keV 
(Fig. 2, top). We used EPIC-MOS data in the range 0.5-2.0 keV 
to verify the presence of a MCD component in the LHS. Fig. 2 
(bottom) shows the data/model ratio (extended to an energy range 
of 5.0 keV) without a blackbody component for the LHS. The best 
fit model requires a disc blackbody with a temperature of ~ 0.22, 
with x^/i" = 1137.8/389. The best fit without a quasi-blackbody 
gives x^ /v = 129617.2/389, which clearly indicates the need for 
this component. Recently, similar results have been obtained for the 
low/hard state of GX 339-4, where an optically thick disc with a 
temperature of ^ 0.2 keV has been reported (Tomsick et al. 2008). 

3.2.2 Simple Model: 0. 7-9. ke V 

We fit both the VHS and LHS simultaneously with a power-law 
plus MCD component. Only the value of A'h was tied between 
the states. The best-fit value found for A^h of 5.170 ± 0.001 x 
10^^ cm~^ is in accordance with that of Kong et al. 2000. We re- 
stricted the value of A'h to 5.1-5.3timesl0^^ for the remainder 
of thus work. Figure 3 shows the spectra with the data/model ra- 
tio. The formally unacceptable fit (x^/u = 10279.4/3860) can be 
seen to be due mainly to the broad iron line and soft energy residu- 
als, and for the VHS, a large Fe Ka absorption edge. 

To provide a physically realistic description of the Fe line re- 
gion we initially added a relativistic Fe line (LAOR, Laor 1991) to 
the MCD and power-law continuum and, for the VHS, a smeared 
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Figure 2. Top: Best fit spectra for the VHS showing the presence of a 
quasi-blackbody with a temperature of Ri 0.76 kcV (see text). Bottom: 
Data/model ratio without a quasi-blackbody component for the LHS. Rev- 
olution 782 and 783 are shown in red and blue respectively (the data are 
combined for plotting purposes only). It is clear from these plots that a 
semi-blackbody component is present in both states. 



edge to phenomenologically model the iron Ka absorption edge. 
The L AOR model describes a broad line from an accretion disc sur- 
rounding a rotating Kerr BH, with an emissivity profile described 
by a power-law of the form ej^) = r^'. The outer disc radius was 
fixed at the maximum allowed value of 400rg. The inner radius of 
the disc, nn, emissivity index, q, disc inclination, i, and the nor- 
malization were free to vary. It should be noted that constraining 
the spin based on the LAOR model, although robust, is only an 
approximation since the identification of nn, as determined from 
LAOR assumes a hard wired spin parameter of a = 0.998. The 
way that the infened BH spin depends on the position of the in- 
ner radius was explored by Dovciak, Karas & Yaqoob (2004), and 
was shown to be consistent with the "true" spin as one consid- 
ers more rapidly rotating black holes (see their Figure 2). We fit 
both the VHS and LHS individually, restricting the value of TVh 
to 5.2 ± 0.1 X 10^^ cm^^. The fit parameters are given in Table 
1. Adding both a LAOR and SMEDGE components significantly 
improved the fit for the VHS, with x^/^ = 2478.1/1652 and the 
LHS with x^/^ = 3069.7/2376 (Fig. 4). For the LHS, residuals 
below 2.0 keV indicates that the simple power-law plus MCD com- 
ponents, predominant in this range, is not a accurate description of 
the continuum, and a more complex reflection component should 
be present. If data below 2.0 keV is removed for this state, and 
the disc temperature, normalisation and column density is frozen. 
Model 1 converges to /'^ ~ 2249.3/2031. The energy range be- 
low 2.0 keV does not affect the direct measurement of the Fe Ka 



Table 1. Results of Fits to XMM-Newton EPIC-pn (VHS) and EPIC-MOS 
(LHS) data. The spectra were analysed in the 0.7-9.0 kcV band. 

Model 1 ( PHABS X SMEDGE X [ PL + DISKBB + LAOR ]) 



Parameter 



Very High State Low/Hard state 



Nu (lO^i cm-2) 

r 

<;r(keV) 

Rmcd 

QLaor 

nn(rg) 
i(deg) 

iJLaor(xlO-3) 
-Esmodge(fcey) 
^smedge 
^^smedge 



5.300-0.004 
2.567 



+0.004 

0.01 
+0.02 
0.01 



2.87: 
0.721 ± 0.001 

2890l 



-10 

-30 



6.97-0.01 



6.82 



+0.03 



i.9i: 

18.2^ 



-0.04 
1+0.02 
-0.01 
,+0.3 
-0.5 

130 ±2 
7.10+0-01 

9 o+O.l 

4.1 ±0.2 



5.100+0006 
1.700-0.002 
0.231 ±0.005 
0.235 ± 0.02 
8100 ± 300 

"•^'-0.06 

O QQ+0.05 

''•^•-'-0.04 
2.8±0.1 

10.0+20 
3.75 ±0.15 



2478.1/1652 3069.7/2376 



line profile since the LAOR function used to model the line only 
extends down to an energy of ~ 3.5 keV (see Fig. 4). We note that 
the parameters values found here differ slightly from those of sim- 
ilar models in Ml and M2 likely due to the restriction imposed on 
the neutral hydrogen column density, Nu and on improved calibra- 
tion. 



3.2.3 More Complex Models: Ver)' High State 

In all our previous fits, the presence of a broad iron emission line 
has been determined and successfully modelled by the LAOR ker- 
nel. The presence of a possible edge at ~ 7.1 keV found for the 
VHS is consistent with that predicted from absorption due to Fe 
A'-shell transition in partially ionised, "warm", material (Ross & 
Fabian 1993; Ross et al. 1996). To date, BHB spectra have been 
modelled by a combination of a line model such as LAOR or 
KERRDISK (Brenneman & Reynolds 2006), a separate reflec- 
tion function such as PEXRIV and a multicolour disc, since no 
self-consistent reflection model had been available. Here, we use 
the reflection model developed by Ross & Fabian (2007, REFHI- 
DEN) to model those components jointly for the VHS. The pa- 
rameters of the model are the number density of hydrogen in the 
illuminated surface layer, H^cn, the value of kT for the blackbody 
entering the surface layer from below, the power-law photon in- 
dex, and the ratio of the total flux illuminating the disc to the to- 
tal blackbody flux emitted by the disc. The disc reflection spec- 
tra is convolved with relativistic blurring kernel KDBLUR, which 
is derived from the code by Laor (1991). The power law index 
of REFHIDEN is tied to that of the hard component. We con- 
strain the value of A^h, inclination, and power-law index to 5.1- 
5.3x10^1 cm"2, 10-20 degrees and 2.5-2.7 as found from Model 
1 above. Using this model we obtained a much improved fit for 
the VHS, with xV^^ = 2348.8/1655 (Model 2). Allowing for 
a broken power-law for the emissivity further improves the fit, 
with xVi' = 2237.8/1653 (Model 3, see Fig. 5 top) and an F- 
Statistic value over the previous fit of 41, implying a probability of 
< 10~i* of a random occurrence. The parameters found for these 
models are shown in Table 2. A direct inspection of Fig. 5 (Top) 
shows the presence of a possible photo-ionisation edge for O VIII 
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Figure 3. Top: Data/model ratio for a simple model consisting of a power-law and MCD components only. The XMM-Newton EPIC-pn data for the VHS is 
shown in black. Spectra for the combined EPIC-MOS (LHS) revolution 782 and 783 are shown in red and blue respectively (the data are combined for plotting 
purposes only). Bottom: Blowup of the EPIC-pn (black) and EPIC-MOS 1 (revolution 782, red) spectrum showing the broad iron line and Fe Ka edge region. 
The data have been re-binned for visual clarity. The EPIC-pn spectrum has been extended to 10 kcV for illustration purposes only. 
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Figure 4. Left: LHS data/model ratio for a model consisting of a simple power-law and MCD components, as well as a LAOR line (Model 1). Spectra for the 
combined EPIC-MOS revolution 782 and 783 are shown in red and blue respectively. The data are combined and re-binned for plotting purposes only. Right: 
Model components for MOS 1 revolution 782. 
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Table 2. Result of more complex fits to XMM-Newton EPIC-pn data for 
GX 339-4 in the Very High State. 



Parameter 


Model 2 


Model 3 


N-a (10^^ 


5.16 ±0.02 


5.17 ±0.03 


r 


2.7-0.01 


2.7-0.02 




1 45+0-25 




kT ( keV) 


0.519tooo6 


0.554 ± 0.004 


//don(xl02lHcm-3) 


4.05 ±0.2 




Illum/BB 


4 4+0-4 


2.0 ±0.1 


^REFHIDEN 


4.410;? 


9 7+0.1 
^■'-0.3 


Qin 


6.84l°-i 


7 6+0-3 
'■"-o.e 


gout 




o 7+0.3 
■^•'-0.8 


^brcak i^g ) 




4 g+0.6 


rin(rg) 


1 sriA+O'OS 
i.B04_Q (JP4 


r, r,q + 0.025 

^•"'-'-0.035 


i (deg) 


19.510-5 


19.9810-02 




2348.8/1655 


2237.8/1653 



Notes.- Model 2 is described in XSPEC as (PHABS X ( KD- 
BLUR X [PL + REFHIDEN ])). Model 3 assumes the accretion disc has a 
broken-power law emissivity profile described by the function KDBLUR2 
in XSPEC. The value of Nn, inclination, and power-law index was con- 
strained to 5.1— 5.3x10^^ cm~2^ 10-20 degrees and 2.5-2.7, respectively 
in both models. 

at ~ 0.86 keV and a narrow line at ~ 6.4 keV. Adding an edge at 
0.86 keV with an optical depth r = S.Slg g x 10"^, as well as a 
narrow Gaussian at 6.4 keV results in x^/*^ = 2144.9/1651 (Fig. 
5 bottom). The presence of the narrow Gaussian component is re- 
quired at the 98 per cent level (F-test probability of 0.02) and can be 
attributed as due to reflection from distant materials. By looking at 
the contribution to (not shown in Figure 5) it is clear that the ma- 
jority of the contribution comes from the energy range 1.6-3.3 keV 
and is likely the result of Au M-shell edges and Si features in the 
detector. Due to the high signal-to-noise ratio achieved, these fea- 
tures are clearly revealed and thus affect the overall x^ fit statistics 
but not the values of the best-fit parameters. 

From the irradiating source luminosity, Lx (defined here be- 
tween 0.2-10.0 keV), disc hydrogen-density (Hdcn) and radius, R 
(approximated to 2.0rin), it is possible to estimate the ionisation pa- 
rameter, (^ = L^/HdenR^), of the VHS as log^ « 4.2, which is 
consistent with previous values found for GX 339-4 in its very high 
state (Ml). Using a self-consistent reflection model we were able 
to constrain the value of the innermost radius to nn= 2.03lQ o35rg. 
This strong constraint can be better appreciated in Fig. 6, where the 
90 per cent confidence lever for nn is shown as the dashed-line in 
the x^ plot, obtained with the "steppar" command in XSPEC. As- 
suming that r,n = Tins we get a spin parameter a/M— 0.9391q qq3 
for the VHS. 

It should be noted that the measured values for the disc tem- 
perature found by the REFHIDEN model is at least 0.15 keV 
less than the value obtained via fits with the additive compo- 
nents (Model 1). This behaviour is as expected (see Merloni, 
Fabian & Ross 2000; Ross & Fabian 2007) since in Model 1 
the soft disc component is modelled with a multi-colour disc 
quasi-blackbody (DISKBB) having a Comptonised surface colour- 
temperature (Tcoi) intrinsically higher than the effective mid-plane 
temperature (T^f /) used in REFHIDEN. The spectral hardening 
factor f^oi = Tcoi/Tefs (Shimura & Takahara 1995) of fa 1.3 
found for the VHS is very similar to the conventional value of 
1.7 ± 0.2 (Shimura & Takahara 1995). It is expected that a sim- 
ilar reduction in the value of kT should occur for the LHS where 




Energy (keV) 

Figure 5. Data/model ratio for the VHS. Top: Model assumes a bro- 
ken power-law emissivity profile and constitute of a power-law and the 
disc reflection function REFHIDEN. The presence of a O Vlll edge at 
0.86 kcVand a possible narrow emission line at » 6.4 kcV can be seen. 
Bottom: Same as above but with an additional narrow Gaussian line and a 
O Vlll edge at 0.86 kcV. The data have been re-binned for visual clarity. 



the effective disc temperature would be less than the value found in 
section 3.3.1 of 0.22 keV. 



3.2.4 More Complex Models: Low/Hard State 

The effective temperature of around 0.15 keV and hydrogen num- 
ber density //den > 2 x lO^^Hcm"^ expected for the LHS fafls 
outside the parameter range of the REFHIDEN model. However, 
the model is being developed and an analyses of the LHS with 
REFHIDEN is left for future work. The disc-blackbody in the LHS 
has a negligible effect on the iron- /fa features above 2.0 keV(see 
Figure 4). For comparison with our results for the VHS, we use RE- 
FLIONX to analyse XMM-Newton EPIC-MOS data for the LHS 
in the range 2.0-10.0 keV with a model similar to REFHIDEN 
but lacking the intrinsic blackbody disc component. The model 
REFLIONX is a revised version of REFLION (Ross & Fabian 
2005) used to describe reflection from accretion disc in AGN sys- 
tems where the blackbody emission is at too low an energy to af- 
fect the Fe Ka emission. It should be stressed that the reflection 
features above 2keV are unlikely to be significantly affected by 
the change from REFHIDEN to REFLIONX. The parameters of 
the model are the iron abundance (set to solar), photon index of 
the illuminating power-law, ionisation parameter, ^, and the nor- 
malization. The disc reflection spectra is convolved with the rela- 
tivistic blurring kernel, KDBLUR. The power law index in RE- 
FLIONX is tied to that of the hard component. We constrain the 
value of the inclination, and power-law index to 10-20 degrees and 
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Figure 6. 'vs plot for GX 339-4 in its Very Higli State. A value 
of i'in= 2.03^Q g3grg is found at the 90 per cent confidence level for 
one parameter, Ax'^ = 2.71 criterion shown by the dashed line. The inset 
shows a close up of the region around the minima. 
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Figure 8. x^ vs rin plot for GX 339-4 in its Low/Hard State. An inner 
radius, ri^= 2.08^q ^JpTg is found at the 90 per cent confidence level for 
one pai'ameter, Ax^ = 2.71 criterion shown by the dashed line. The inset 
shows how behaves closer to the minima. The inner radius found for the 
VHS is shown in red for comparison. 



2; _ 




Energy (keV) 

Figure 7. LHS data/model ratio for a relativistic blurred disc reflection 
model (REFLIONX) and power-law as described in the text. The spectra 
for the combined £P/C-MOS revolution 782 and 783 are shown in red and 
blue respectively. The data have been combined and re-binned for plotting 
purposes only. 



1.4-1.6, respectively. The hydrogen column density, A'h is fixed at 
5.17 X 10^^ cm^^, the best fit value found for the VHS, as we do 
not expect it to vary. The best fit obtained with the blurred RE- 
FLIONX model is shown in Fig. 7 and detailed in Table 3. This 
model gives jv = 2242.5/2031 with an emissivity index of 
3.065 ± 0.05 indicating a standard "lamp-post" emissivity pro- 
file. The value found for the inclination of 2O-1.7 deg. is in agree- 
ment with that for the VHS. The low disc ionisation parameter, 
log(5)R; 3, is consistent with that expected for low disc temper- 
atures. At 90 per cent confidence, this model gives constraint on 



the innermost stable radius of rin= 2.1 



(see Fig. 8). If we 



include the energy band 0.7-2.0 keV to the above fit, a large low en- 
ergy residual is present due to the disc emission. By modeling this 
with a DISKBB component, a best fit of x^ jv = 3070.9/2388 is 
achieved in the full 0.7-10.0 keV range with a disc temperature of 
0.201 ± 0.003 keV as in 53.2.1. 



Table 3. Results of fits to XMM-Newton EPIC-MOS data for the Low/Hard 
state of GX 339-4. The spectra was analysed in the 2.0-10.0 kcV band. A^h 
was fixed at the quoted value. 



Parameter 



Low/Hard State 



Nu (10^ 



5.17 



r 
1 

nn(rg) 
i (deg) 

5(erg cms~^) 
^Jreflionx(xIO-S) 



1 48+°-°^ 
0.11 ±0.01 
3.065 ±0.05 
2.08l»:lS 

20-1.7 
1350 ± 100 

4.0j 
2242.5/2031 



-1+0.4 
-1.0 



3.3 Joint XMM-Newton - RXTE spectrum analysis 
3.3.1 Very High State: 0.7-lOO.OkeV 

In order to check the robustness of our results, we extended the 
fit from the EPIC-pn spectrum to include the energy range 0.7- 
100 keV, using RXTE data. PCU-2 data was fitted between 8.0- 
25.0 keV. This resulted in a poor fit, with x^/f = 4017.0/1727. 
It should be stressed that the quality of the EPIC-pn data far out- 
weights that of RXTE and thus a statistically worst fit is inevitable 
in the full range. However, most of the residuals are accounted for 
by allowing the power-law index to vary. The best fit value of F = 
2.583 ± 0.007 is in accordance to that found in section 3.1 for the 
RXTE continuum. The parameters for fits to the combined EPIC- 
pn and RXTE spectrum are listed in Table 4 and shown in Fig. 9 
(left). It is clear that the model is a very good description of the 
spectrum (see Fig. 9), with x^ /i' = 2549.3.0/1718 in the full 0.7- 
100.0 keV range. Most importantly, the value for the inner radius. 



Tin found here of 2.02 
3.2.3. 



+0.02 



is similar to that found in section 
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Figure 9. XMM-Newton and RXTE spectra of GX 339-4 fit jointly with a disc reflection model convolved with the relativistic blurring kernel, KDBLUR. 
Left: VHS spectrum in the 0.7-100.0 keV range. Right: LHS spectrum in the 2.0-100.0 keV range. The best fit indicates a value for r-^^ of 2.02l^Q yg rg and 
2.04^0 QjTg for the VHS and LHS respectively. £P/C-MOS/pn data are shown in black. RXTEVCV-l and HEXTE daXa are shown in red and blue respectively. 
The data have been re-binned for visual clarity. 



Table 4. Joint XMM-Newton and RXTE spectral fits with relativistic bluired 
disc reflection models. 



Parameter 


Very High State 


Low/Hard State 




5^^00+0.004 


5.17 


r 


2.583 ± 0.007 


1 ..,+0.005 
-^•"*'^-0.02 


Rpi, 


^•"^-0.09 


0q4+0'005 


kT{ keV) 


0.585 ±0.001 




/fden(xl02lHcm-3) 


6.6 ±0.2 




Illum/BB 


1.00 ±0.02 




^(erg cms~^) 


> 10000 


13301^0 


■Rrefhiden 


1 02+0 02 
^■"^-0.06 


Rreflionx{10^^) 




4.4 ±0.2 


Qin 


7 05+0-05 


3.16 ±0.05 


9out 


3.0+01 




^brcak i^'g ) 


6.01°:? 




nnCr-g) 


9 r,9+0.02 

^■'^^-o.oe 


^■^^-0.02 


(■ (deg) 


20.0-0.3 


20.0-1.3 




2549.3/1718 


2316.6/2095 



Notes.- VHS was modelled with REFHIDEN and a broken power-law 
emissivity profile, as described in the text. For the LHS, the disc reflec- 
tion model REFLIONX was used, and the spectra was fitted in the range 
2-100 keV. 

3. 3. 2 Low/Hard State: 2-1 00 ke V 

A similar extension on the energy range of the LHS was made, 
with RXTE data being used in conjunction with EPIC-MOS. A best 
fit of x^/j^ = 2316.6/2095 over the full 2.0-100.0 keV energy 
band is found. The various parameters are shown in Table 4 and the 
resulting spectra in Fig. 9 (right). The value for the inner radius, nn 



found here of 2.04_ 



4 DISCUSSION 



, Tg is similar to that found in section 3.2.4. 



The spectral modelling of the VHS suggests that the surface layer 
of the accretion disc is highly ionised, with ^ ~ 10'* erg cm s^^. 
In this state, the iron in the outer layer of the disc is fully ionised 
and regions > 2 Thomson mean free paths below the surface pro- 
duces a strong iron-K absorption edge. Narrow Ka line emission 



from this region is then Compton-broadened as it scatters out of 
the disc. The strong presence of the iron-K edge in the VHS can 
be seen in Fig. 3 and quantitatively appreciated by the high op- 
tical depth (r — 2.3to g) '^^ (phenomenological) component 
SMEDGE in Model 1. The best-fit REFHIDEN model clearly 
shows the large K-shell absorption feature and weak Ka emis- 
sion line characteristic of the VHS in GX 339-4 (Fig. 10, left). 
Note that similar features have also been observed in the VHS of 
Cygnus X-1 (Done et al. 1992). The ionisation parameter found for 
the LHS, ^ — 1350 ± 100 erg cm s^^, is consistent with the disc 
being moderately ionised and having a low apparent temperature. 
In this state, the illuminated accretion disc results in a strong Fe- 
Ka line emission from the top layers. Compton-broadening, al- 
though present, cannot explain the highly broadened and skewed 
line shape (see Fig. 3), where the low energy red wing extend 
down to f» 4 keV. Figure 10 (right) shows the best-fit REFLIONX 
model, prior to (top) and after blurring (bottom), for the LHS. As 
opposed to the VHS, in the low/hard state the Fe-Ka line is clearly 
seen. The value for the ionisation parameter for the low/hard state 
of 5 ~ 10* ergcms"^ reported by Tomsick et al. (2008) is an or- 
der of magnitude higher than the present result. At these values, 
the iron is fully ionised and should not produce an iron Ka reflec- 
tion line (Matt, Fabian & Ross 1993; Young, Ross & Fabian 1998). 
The apparent inconsistency in their results can be attributed to the 
use of PEXRIV (Magdziarz & Zdziarski 1995) as the reflection 
model. This model does not account for diffusion of photons in the 
disc and thus rectify broadening caused by Comptonisation by in- 
creasing the ionisation parameter. 

The obvious differences in the resulting spectra of the two 
states can be ascribed to the different ionisation states of the disc. 
Previous attempts to model the spectra of Galactic BHB used the 
relativistic blurring of the Ka line to obtain the innermost radius, 
Tin and thus the spin parameter. In the present work, the full re- 
flection spectra for the two extreme states was convolved, and the 
bluiTing parameters were obtained not just from the Ka line but 
from all of the reflection features. This is particularly important in 
the case of the VHS, where the ¥e-Ka emission line is not the 
dominant feature of the reflection component and Compton scat- 
tering needs to be fully accounted in the reflection model. In this 
state, a steep inner disc emissivity index of gin = 7.6lQ g, within a 
radius rbroak = 4.9lo ^'"g required, indicating that the corona is 
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Energy (keV) Energy (keV) 

Figure 10. Top: Best-fit model components prior to blurring for the VHS {left) and LHS {right). The combined model, reflection, and power-law components 
are shown in black, red and blue respectively. Iron is fully ionised in the VHS, thus resulting in a very small ¥e-Ka line emission and a large iron-K absolution 
edge. In the LHS, moderate ionisation means that the contribution to the spectra from Fe-Zfa hne emission is greater (see text). For the LHS, we only show 
EPIC-MOSl revolution 782 for visual clarity. Bottom: Same as above, with relativistic blurring. 



centrally concentrated. The model constrains the inner radius of the 
accretion disc to rin= 2.03to Q35 7g at the 90 per cent confidence 
level. It should be noted that the value for rin quoted above for the 
XMM-Newton observation is consistent with that found for the full 
full XMM-Newton plus RXTE fits, indicating that the model is an 
accurate description of both the reflection features as well as the un- 
derlying continuum. Assuming that emission within the innermost 
stable orbit is negligible (see Reynolds & Fabian 2007), the value 
of Tin found here translates to a black hole spin of 0.939lg'Q03 for 
the VHS. In the LHS, spectral fitting using the model REFLIONX 
resulted in an inner radius of rin= 2.08l:J5 ioTg. This translates to 
a spin parameter of O.93t5!:o2^ for the LHS. 

The value for the spin parameter found for both states of 
GX 339-4 are within one per cent of one another and falls within 
one sigma error. It has been argued that bleeding of the iron 
line emission to regions inside the innermost stable radius may 
cause systematic errors in the derived value of the spin param- 
eter (Reynolds & Begelman 1997; Krolik 1999). Using a high- 
resolution 3-D MHD simulation of a geometrically-thin accretion 
disc, Reynolds & Fabian (2007) have shown that the ionisation 
edge is within ~ O.Sr-g of the innermost stable circular orbit for 
a non-spinning BH. However, it was shown by the same authors 
that this bleeding decreases as the position of the innermost radius 
approaches the event horizon, and hence the spin inferred from the 
position of rin becomes much closer to the true spin as one consid- 
ers more rapidly rotating black holes. Similar results were reported 



by Dovciak, Karas & Yaqoob (2004, see their Figure 2). In order 
to verify our results against any systematic variation, we modelled 
the VHS with a KERRDISK line profile (Brenneman & Reynolds 
2006). The spin, which is a free parameter of the model, was found 
to be 0.952t°;[J°^, lying within ~ 1 per cent of the value inferred 
from the innermost radius. 

The reflection model, REFHIDEN, assumes a single- 
temperature accretion disc. Although this is not a realistic claim, 
we believe it unlikely to have any significant effect on the inferred 
innermost radius of emission. In order to verify this hypothesis, we 
approximated a "real" disc with inner radius increasing logarithmi- 
cally from 2rg to 6.78rg. Each point on the disc was assumed to 
radiate like a blackbody with an effective temperature that scales 
with radius as r~"^^* starting at 0.9 keV. Within this region, the il- 
luminating flux scaled as r~®. Using the EPIC-pn response file, we 
modelled 1 ks of simulated data with a single temperature REFHI- 
DEN. As expected, the model constrained the various parameter, 
with an effective temperature of ~ 0.52 keV and an inner radius 
of 2.041^Q Q2o ''g- ^ further check on any inconsistency that 
may arise from using a single temperature disc reflection model 
to constrain the spin of the black hole, we investigated the VHS 
with a different thermal model, (KERRBB, Li et al. 2005), which 
includes relativistic smearing in a disc with radial temperature gra- 
dient. The black hole spin, a free parameter in the KERRBB model, 
was found to be 0.93 ± 0.02, consistent with that inferred from the 
single temperature REFHIDEN model. 



10 R. C. Reis et al. 



5 CONCLUSIONS 

We have analysed XMM-Newton spectra of GX 339-4 in both its 
very high and low/hard states. Looking at the difference in the spin 
parameter between the two states, as well as that derived from the 
various independent models for the VHS, we can estimate the sys- 
tematic error in the iron line method to be about 1 per cent. By 
using a reflection model which intrinsically accounts for Compton- 
isation and blackbody emission, we infer that the spin parameter of 
GX 339-4 is 0.935 ± 0.01 (statistical) ± 0.01 (systematic). 
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